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Abstract 

It is shown that spatially homogeneous solutions of the Einstein 
equations coupled to a nonlinear scalar field and other matter exhibit 
accelerated expansion at late times for a wide variety of potentials 
V. These potentials are strictly positive but tend to zero at infinity. 
They satisfy restrictions on V /V and V" /V related to the slow-roll 
approximation. These results generalize Wald's theorem for space- 
times with positive cosmological constant to those with accelerated 
expansion driven by potentials belonging to a large class. 

1 Introduction 

In recent years astronomical observations have shown the necessity of incor- 
porating accelerated expansion into the description of the evolution of our 
universe. Unfortunately there is no clear understanding of the mechanism 
leading to the acceleration. The simplest possibilities are a cosmological con- 
stant or a nonlinear minimally coupled scalar field with a potential. There 
is now a large and growing literature on nonlinear scalar fields with different 
potentials and many more exotic models. The observational data available 
at present only gives weak criteria for deciding between different models. For 
instance it is not even clear that it can rule out the case of a cosmological 
constant (cf. pQ). It is to be hoped that this situation will improve due the 
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accumulation of more data in the next few years. In the meantime it makes 
sense to advance our theoretical understanding of the models available. 

The present paper is concerned with the dynamics of spatially homoge- 
neous solutions of the Einstein equations in the presence of a minimally cou- 
pled nonlinear scalar field and other matter. The latter represents ordinary 
baryonic matter and non-baryonic dark matter. In the literature on models 
with accelerated cosmological expansion one or more of the following effects 
are often ignored: shear, ordinary matter, spatial curvature. The usual in- 
tuitive picture says that they will be negligible under suitable circumstances 
in a phase of accelerated expansion. A rigorous proof that this is true was 
given by Wald |23| in the case the acceleration is caused by a cosmological 
constant. The first main result of the present paper is that Wald's theorem 
can be generalized to the case where the cosmological constant is replaced 
by a scalar field belonging to a large class. This class apparently includes all 
of the potentials considered in the literature for which accelerated expansion 
is expected to continue indefinitely. 

A concept which is frequently used when studying the dynamics of models 
with accelerated expansion is that of the slow-roll approximation. It says 
that under certain circumstances the equation of motion of the scalar field 
can be replaced by a simpler equation while retaining the essential qualitative 
features of the evolution. The precise meaning of the slow-roll approximation 
has been investigated in [T7]. It was shown that if a solution approaches 
the 'slow-roll attractor' many conclusions follow. The second main result of 
the present paper is a proof that for a wide class of potentials all solutions 
approach the slow-roll attractor. 

Since the potential of the scalar field leading to accelerated expansion is 
not known it is useful to try to analyse the dynamics of spacetimes with a 
nonlinear scalar field for large classes of potentials assumed only to satisfy 
some qualitative restrictions. This programme was carried out in [22] for 
potentials with a strictly positive lower bound and solutions of Bianchi types 
I- VIII. The spacetimes analysed may contain matter satisfying the dominant 
and strong energy conditions in addition to the scalar field and the results 
are a direct generalization of Wald's theorem [33] 011 spacetimes with positive 
cosmological constant. As pointed out in [22] the next natural case to look 
at in order of increasing difficulty is that where the potential is everywhere 
positive but may tend to zero as the scalar field <fi tends to ±oo. Furthermore 
it should be assumed that the fall-off at infinity is not too fast, since otherwise 
accelerated expansion may not occur. 
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To see the need for a restriction on the fall-off of the potential it is useful 
to look at explicit solutions given by Halliwell jT^j. These solutions are 
homogeneous and isotropic and spatially flat. Their only matter content is 
a scalar field with an exponential potential. In notation which is convenient 
for the following the potential is given by V(<p) = Voe~^^ x ^ ', where A is a 
positive constant. If these solutions are written in Gauss coordinates the 
spatial metric is of the form gij(t) = a 2 (t)Sij where the scale factor a(t) is 
proportional to t 2 ^ x . The expansion is accelerated iff the power occurring 
is greater than one, i.e. iff A < \/2. Thus if the aim is to consider models 
with accelerated expansion an assumption should be made on the potential 
which rules out too rapid exponential decay. This can be done by means of 
a restriction on V'/V. 

Acceleration due to an exponential potential is often called power-law 
inflation since the scale factor has a power-law behaviour in the simplest 
cases. A class of models where the scale factor behaves in a way intermediate 
between a power law and the exponential expansion resulting from a positive 
cosmological constant is associated with name intermediate inflation The 
specific case considered in [H] is where the scale factor is proportional to 
exp(t / ) for < / < 1. This can be obtained from a potential which looks 
asymptotically like a negative power but is not exactly a power. Note that 
if a cosmological constant is thought of as a constant potential the potential 
in intermediate inflation is intermediate between the cosmological constant 
and the potential for power-law inflation as far as its asymptotic behaviour 
for large <fi is concerned. The aim of this paper is to analyse the case of 
general potentials which are positive, tend to zero as (f) —>■ oo and are such 
that —V'/V satisfies an upper bound which rules out exponential potentials 
with A > 

Under an additional assumption on V" /V further information is obtained 
on the asymptotic behaviour of the solutions. These properties are related 
to the slow-roll approximation ^7] . It is shown that the results obtained are 
general enough to apply to many potentials considered in the literature. 

After some necessary background material has been presented in section 2, 
the slow-roll approximation is discussed in section 3 for spatially flat isotropic 
models with a non-linear scalar field as the only matter content. Section 4 
deals with the generalization to all models of Bianchi types I- VIII and to the 
case where arbitrary matter fields satisfying the dominant and strong energy 
conditions are present in addition to the scalar field. In section 5 it is shown 
that the results which have been obtained are applicable to many potentials 
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considered in the literature. The late-time dynamics of perfect fluids and 
collisionless matter is treated in section 6. Section 7 sums up what has been 
achieved and indicates some avenues for future progress. 

2 Background 

This section recalls the equations which are needed in the following and some 
basic facts which were proved in [221 • Consider a spacetime with vanishing 
cosmological constant which contains a nonlinear scalar field and some other 
matter. Suppose that the other matter satisfies the dominant and strong 
energy conditions. The energy momentum tensor is 

9a/3- (1) 

where T*| is the energy-momentum tensor of the matter other than the scalar 
field. Assume that the potential V is C 2 and non-negative. Restricting to a 
Bianchi spacetime leads to the following basic equations: 

-4VT0 2 - l -a ab a ab + l -R- A*{p M + UrS M ) (2) 
-3H4> - V'{(j>) (3) 
|[0 2 + 2V&)) + l(a ab a ab -R) + fp M (4) 

Here H is the Hubble parameter, a ab is the tracefree part of the second 
fundamental form, R is the scalar curvature of the spatial metric and p M 
and tr5' M are projections of the energy-momentum tensor T£i. For Bianchi 
types I to VIII the inequality R < holds. Assuming that H > at some 
time (so that the model is expanding at that time) it follows that H > at 
all times. For more details see [2*2*] . 

If H is positive at some time then it is everywhere positive and non- 
increasing and so tends to some limit Hi > 0. In Theorem 1 of section 3 of 
j22] some general results on late-time dynamics for the above equations were 
proved unter the following three assumptions: 

1- V(4>) > Vo for a constant Vo > 

2. V is bounded on any interval on which V is bounded 
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dH 
~dt 

= 

H 2 = 



3. V tends to a limit, finite or infinite, as tends to oo or — oo. 

Suppose for a moment that the first of these assumptions is dropped but the 
other two are retained. The first part of the proof of Theorem 1 of |22] does 
not use assumption 1. and so it can still be applied. The conclusion is that 
if a spacetime of the type described above exists globally in the future then 
0-^0 and V(<f>) tends to a limit V\ > as t — > oo. If Hi > then the rest 
of the proof of the theorem also goes through and the same conclusions are 
obtained as in that theorem. In particular the expansion is accelerated at 
late times, the model isotropizes and the contribution of matter and spatial 
curvature to the field equations becomes negligible. Moreover 3H 2 /8nV — > 1 
and V'(<f>) — > as t — ► oo. 

It remains to consider the case Hi — 0. In that case it follows that —>■ 
and V(<f>) —>■ as t — > oo without having to make the additional assumptions 
2. and 3. When Hi = it follows that a ab a ab , R and p M converge to zero as 
t — > oo. The quantity 

Z = 9H 2 - 24tt[0 2 /2 + V{4>)} = ^{a ab a ab - R) + 2Ati P m (5) 

satifies the inequality dZ/dt < —2HZ. In |22] this was combined with the 
fact that H has a positive lower bound to show that Z decayed exponentially. 
In the present more general case this is no longer true but the differential 
inequality can be used in a different way. In fact H 2 > Z/9 and so dZ/dt < 
— |Z 3 / 2 . It follows that Z = 0(t~ 2 ) and this gives decay rates for a ab a ab , R 
and p M . Note that the conclusions up to this point also hold in the absence 
of a scalar field. Assume now that V is positive. (This assumption will be 
maintained for the rest of the paper.) Then V(<f)(t)) — > and so — ► oo 
or0— > — ooast-^oo. These cases can be transformed into each other by 
reversing the sign of and so it can be assumed without loss of generality 
that — > oo. In this case V(0) — > as — > oo. 

3 The slow-roll approximation 

In this section we specialize to the isotropic and spatially flat case with a 
nonlinear scalar field as the only matter. This is the situation most frequently 
considered in the literature. In this case H — a/ a where a is the scale 
factor. It will be shown that certain assumptions on the potential V imply 
accelerated expansion at late times. The assumptions are 
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1. V(<f>) > with V(4>) -»• as -> oo 

2. < 

3. V r, (0)/^(0) -> as -f oo. 

Consider a situation where, as in the last section, — > oo as t — > oo. 
Note that > for t sufficiently large. For must be positive at some 
time. After that it must remain positive since = and assumption 2. 
above imply that > 0. When it is known that > the assumptions of 
the last section can be weakened while maintaining the same conclusions. 
Assumption 2. only needs to be required for intervals whose left endpoint 
is finite. In assumption 3. only the existence of the limit for t — > oo need 
be required. These weakened versions of assumptions 2. and 3. of the last 
section follow from the assumptions 1.-3. of this section. The following 
argument was inspired by a paper of Muslimov [TH] on the isotropic case but 
we will not use the complicated transformation of variables of that paper 
directly. Since is monotone we can regard functions of t as functions of 0. 
The key relation is the following: 

d (3H 2 \ _ 3H 2 ( 8tt0 V'\ 
d4> \8^VJ ~ 8W \ If ~ V) ^' 

Theorem 1 Consider a spatially flat homogeneous and isotropic solution of 
the Einstein equations coupled to a nonlinear scalar field with potential V 
of class C 2 satisfying conditions 1.-3. above. Suppose that the solution is 
initially expanding (H > 0) and that > at some time. If the solution 
exists globally to the future then 3H 2 /8irV(4>) — > 1 as t —>■ oo and a > for 
t sufficiently large. 

Proof In order to control the first term on the right hand side of (jUJ) note 
that 

_ [3 2V^) 

and that the right hand side is an increasing function of 3H 2 /8itV. Moreover 
it tends to zero as 3H 2 /8ttV — * 1. Let C\ > 8ir/3 be a constant and let 

C 2 = 8-k^3/ 'An — 2/C±. Since V'/V — > as — > oo there exists some 0i 
such that — V'(0)/V(0) < C2/2 for all > 0i. At any point where > 0i 
and H 2 /V > C\ the expression in brackets on the right hand side of (jBJ) is 
bounded above by —C2/2. Hence H 2 /V, if it is ever greater than Ci, reaches 



6 



C\ at some greater value 02 of and remains below that value for all > 02- 
Because of the freedom to choose Ci, and the fact that 3H 2 /8ttV > 1, it 
follows that 3H 2 /8nV — > 1 as —>■ oo and hence as t —>■ oo. Thus the first 
conclusion of the theorem has been proved. 

The inequality H + H 2 > is equivalent to a > and is thus the criterion 
for accelerated expansion. Substituting for H and H using the field equations 
shows that it is equivalent to the condition that 3H 2 /8nV((f)) < 3/2, which 
is fulfilled at late times, giving the remaining conclusion of the theorem. 

Remark The inequality H + H 2 > implies that H(t) > (t + C)^ 1 for some 
C > and t large and that J t ^° H{t)dt = oo for any t . 

An analogous result can be obtained in the case that assumption 3. above 
is replaced by the weaker condition that limsup(— V'/V) < a for a suitable 
constant a. 

Theorem 2 Consider a spatially flat homogeneous and isotropic solution of 
the Einstein equations coupled to a nonlinear scalar field with potential V of 
class C 2 satisfying conditions 1. - 2. above with a = limsup(— V'/V) < 4^/7r. 
Suppose that the solution is initially expanding (H > 0) and that > at 
some time. If the solution exists globally to the future then a > for t 
sufficiently large. 

Proof Suppose that C\ > (3 where 



(3 4vr 64vr 2 v ' 

and define Ci in terms of C\ as above. Then C2 > 01. There exists some 
0! such that -V'/V < (C 2 + a)/2 for > fa. If > X and H 2 /V > C x 
then Sinp/H — V'/V < (a — C 2 )/2 < 0. It follows that, as in the previous 
case, (H((f))) 2 /V(4>), if it is ever greater than Ci, reaches C\ at some greater 
value 02 of and remains below that value for all > 02- Hence in the 
limit — > 00 we have limsup3if 2 /87rl / < 3/3 /8tt. Thus the condition for 
accelerated expansion is 3f3/8n < 3/2, i.e. a < iy/n. In particular the 
inequality for accelerated expansion in the case of an exponential potential 
with A = a/y/8n is recovered. 

Remark For an exponential potential it follows from ^Sj that 3H 2 /8ttV — > 
3f3/8n. 

In the case that V'/V — > as t — * 00 it is possible to obtain further in- 
formation on the asymptotic behaviour of the solutions which relates to the 
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so-called slow-roll approximation ^7]. This says that under the assump- 
tion on the potential that V" /V is bounded the term in the equation 
of motion of the scalar field becomes negligible at late times. Equivalently 
3H<p/V — > — 1. In the terminology of this shows that for a potential 
satisfying this assumption all solutions approach the slow-roll attractor. 
Theorem 3 Let the assumptions of Theorem 1 hold and assume in addition 
that V"/V is bounded as <p -> oo. Then 3H<p/V -> -1 and 4>/3H<j) -> as 
t — >■ oo. 

Proof The desired limiting behaviour can be proved using the relation 




-H 



3H6\ H6 

1 + — ) + v^ 




(9) 



Let C\ > be a constant. There exists a time t\ such that for t > t x the 
inequality (j)/H)(V"/V) > -C x /(1 + C x ) holds. If 3#0/V < -1 - C x at 
some time t > t x then it follows from Q that 




2C X H fHf 



(i + Ci) V v 



(10) 



Since if is not integrable it follows that 1 + 3H(p/V' must reach —C x after 
finite time and cannot become less than this value again. Using the freedom 
to choose C x shows that liminf(l + H<f>/V) > 0. In particular, since V is 
negative, it follows that H<f>/V is bounded. Now choose t 2 such that the 
modulus of the second term in the square bracket in ((OJ) is less than C x /2 for 
t > tz- Suppose that at some time t > t% the inequality 1 + 3H<p/V' > C x 
holds. Then 1 + 3H(p/V is decreasing there. Moreover the rate of decrease 
is such it must reach C x in finite time. It follows that 1 + 3H<p/V — > as 
t — > oo. The equation of motion for <fr implies that <p/3H(j) — > 0. 

Remark For an exponential potential it follows from ^H] that 3H<p/V — > 
-(1 - A76)- 1 as t -> oo. 



4 The anisotropic case with matter 

Here we generalize the results of the last section to anisotropic solutions with 
matter. 

Theorem 4 Consider a solution of the Einstein equations of Bianchi type I- 
VIII coupled to a nonlinear scalar field with potential V of class C 2 satisfying 
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conditions 1. - 3. of section 3 and other matter satisfying the dominant and 
strong energy conditions. Suppose that the solution is initially expanding 
(H > 0) and that > at some time. Then 

a) if the solution exists globally to the future then 3H 2 /8irV(<f)) — ► 1, cx a b<J ab /H 2 , 
R/H 2 , and p M / H 2 tend to zero as t — > oo and H + H 2 > for t sufficiently 
large and 

b) if in addition V" /V is bounded then 3H<p/V — > — 1 as t — > oo. 
Proof Equation © is replaced by 



d (3H 2 



8irV 



3tt0 V 

IT ~ 17 

#0 V2 6 V 3 



3H 2 
8ttV 
2 



while (|ZJ) is replaced by 



(11) 



H 



3 

Air 



2V(<f)) 1 a ah a ab -R 2p 



Ml 1/2 



if 2 



57T 



# 2 



# 2 



(12) 



It will be shown that (jllj) implies that 0/i? — * as — > oo. Given a 
constant Ci > there exists 0i such that V'/V < C\/2 for all > 0i. Hence 
if j>/H > Ci/8tt and > 0i then d/d(j)(3H 2 /8nV) < -{C X /2){3H 2 /&<kV). 
Thus eventually 0/i? reaches the value C\ and does not return. This proves 
the desired result. Now 



H 



Hep \2 



1 



ah 



-R + Atx[ P 



< -V487T + VA8n 



R 
H 2 \ 6 



V48 



71 



3 \H 



M 



{ 3H 2 J 



H 2 



lib 



in 



f' sM 

1 



o 



3H 2 



v / 48Tp M 



-i 



H 



(13) 



3p M + trS' M 



Provided (j>/H < 1/Vl2n then 
d (3H 2 



8nV 



< 



V48n 1 



3H 2 



V 



'3H 2 ' 
8nV 



(14) 
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Thus arguing as in the isotropic case shows that 3H 2 /8nV — > 1 as — > oo, 
i.e. as t — > oo. It follows that a ab a ab /H 2 , R/H 2 and p M /H 2 tend to zero 
as t — > oo. In particular if + H 2 > at late times and there is accelerated 
expansion. 

The argument to obtain the slow-roll approximation works just as in the 
isotropic case. Equation © is changed by replacing 4n<j) 2 in the last bracket 
by 



47T0 2 + \a ah a ab -U + An (p M + \trS M 
I o V o 



(15) 



The extra terms which are added have the right sign and decay as t — > oo 
and this is all that is needed. 

With some more work this result can be strengthened. Let e > 0. It 
follows from the theorem that 

H + eH 2 = H 2 (e-4ir0 2 /H 2 + o(l)) (16) 

Using what is known about (j)/H shows that for t large the right hand side 
is positive. Thus a differential inequality is obtained which implies that 
H > e (C + t)^ 1 for a constant C > 0. Thus H falls off slower than any 
positive multiple of t~ x . It follows from the inequality dZ/dt < —2HZ that 
Z decays faster than any power of t. Hence Z/H 2 decays faster than any 
power of t. In particular a ab a ab /H 2 decays faster than any power of t. 



5 Applications to potentials in the literature 

In the literature on inflation and quintessence many choices of potential have 
been considered and we will not attempt to give a comprehensive survey. 
Large classes of potentials where the theorems of this paper are applicable 
will be identified and a number of examples in the literature belonging to 
these classes will be pointed out. Then it will be shown by example how the 
theorems can be used to prove asymptotic expansions for the solutions. 
Consider the following class of potentials: 

V{<t>) = U o (log0)^exp(-A0 m ) (17) 

The case p = was studied in while the cases m = and n = 0, m = 1 
were studied in The requirement that V — > as <fi — > oo leads to some 
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restrictions. Suppose first that A < 0. Then it must be assumed that m < 0. 
If m = it must further be assumed that n < and if n = it must be 
assumed that p < 0. If A > it must be assumed that m > or that n and p 
satisfy the restrictions already listed. If A = then n and p must also satisfy 
those restrictions. To compare with the other hypotheses of the theorems 
the first and second derivatives of V must be computed. 



H0) 

V"{<j>) 



p 



'log. 



+ 



n 



Xmcf) 



m—l 



V 



is) 



p(p — 1) p(2n — 1) n(n — 1) 



-Am 



2 (log0) 2 2 log0 
2p m + 2n — 1 



>log( 



+ 



y (19) 



Evidently V'/V tends to zero as — > oo iff m < 1. When m—l the quantity 
—V'/V converges to a positive limit. It can easily be checked that under the 
restrictions already made V is negative for <p sufficiently large. In the case 
that < m < 1 and A > it follows that V" /V behaves asymptotically like 
— Am0 m_1 and so Theorem 4 is applicable. In the case where m < 0, n < 
the quantity V" /V behaves asymptotically like (n — l)/0 and so Theorem 
4 again applies. The case where m < 0, n = and p < is similar, since 
there V"/V behaves asymptotically like (p — l)/01og0. When m = 1 the 
assumptions of Theorem 4 are not satisfied but it follows from Theorem 2 
that in the spatially flat isotropic case with scalar field alone the models 
exhibit accelerated expansion. 

Suppose now that W is a potential satisfying some of the assumptions 
of the theorems and let V be a potential for which V = W(l + o(l)), 
V' = W'(l+o(l)) and V" = W"(l+o(l)) as <p -> oo. Then V satisfies the cor- 
responding assumptions. In particular limsup(— V'/V) = limsup(— W'/W). 
As an example, consider the potentials for intermediate inflation given in [H] . 
These are of the form: 



V{<j>) = V o (0 - 0o) n + W - 0o) n ~ 2 (20) 

for constants n, <po, Vq and V\ with n negative and Vo positive. Taking 
W = Vo0 ra allows the above observation to be applied in this case to see 
that Theorem 4 applies. Another example where this procedure works is in 
the case V = V (exp M p /(j> - 1) [Hj. Take W = V M p /<f>. The case where 
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limsup(— V'/V) is positive gives a proof of late-time accelerated expansion 
for spatially flat isotropic models with scalar field alone in a number of cases. 
There follow some examples taken from a list of potentials in [21] . If V = 
Vb(coshA0 - l) p E3 with -Ay/n < pX < then take W = (V /2P)e pX ^. If 
V = V r o sinh" Q (A0) [22] with < aX < A^ take W = {2 a V )e- aX<l} '. If V = 
Voie^ + e 13 ^) with -Ay/n < an < and (3k < an then take W = V e aK ^. 
If V = V [(<p-B) a + A]e- X<t> and A > then take W = V (p a e- X <t> for a > 
and W = 1/o-4e- A< ^ for a < 0. 

In asymptotic expansions were written down for solutions correspond- 
ing to certain potentials. With the theorems of this paper in hand these 
asymptotic expansions can be given an interpretation which can be proved 
rigorously. For simplicity only the case p = 0,n = 0,0<m<l,A = l 
of (fTTj) will be discussed. Under the hypotheses of Theorem 4, with this 
choice of potential, 3H 2 /8nV — > 1 and 3H<p/V — > — 1 as t — » oo. It can be 
concluded that the solutions satisfy 



as t — ► oo. This should be compared with equation (2.9) of [20]. Substituting 
the specific form of the potential gives 




(21) 




(22) 



This can integrated using the identity 




(23) 



to give 




(24) 



This implies that 




(25) 



Taking logarithms gives 



\og4> 



(log log t + log 2) + o(l) 



(26) 
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Substituting this back into the expression for m gives 

4(1 — m) . . / m 2 I Vr 



21ogt- 



rn 



- (log log t + log 2) + 2 log 



o 



2 V 24vr 



+ o(l) (27) 



Now the asymptotic form of the potential can be computed. 

2 



V{ct>{t)) = Vo 



log 



m 2 / Vr 



o 



2 V 24vr 



r 2 (21ogt) 4(1 - m)/m (l + o(l)) (28) 



From this the asymptotic form of H can immediately be computed H. 



i/(/) = (S^r„/3) 1 / 2 to g [^J^|r 1 (2i«,^) 2,1 -" ,,> (i4 -(D) (201 



To end the section an example will be presented of a potential which does 
not satisfy the assumptions of Theorem 4. This is given by: 



1/(0) = 0" 1 + gin (02) 



(30) 



In this case V'{4>) < for large, V'/V -> as t -> oo, V"/V -> as t -> oo 
but V" /V is unbounded as t — > oo. 



6 Specific matter models - perfect fluids and 
collisionless matter 

All the theorems of the previous sections have been obtained under the as- 
sumption that a solution exists globally in the future. If a suitable well- 
behaved model is chosen for the matter other than the scalar field then global 
existence for arbitrary initial data can be proved. In the case of a perfect 
fluid with a linear equation of state this has been done in [22] • In the case 
of collisionless matter described by the Vlasov equation it has been done in 

Once global existence has been shown it is possible to obtain detailed 
information about the asymptotics of the matter fields. Consider an untilted 
perfect fluid with linear equation of state p = (7 — l)p. The state of the 
fluid is described entirely by its energy density [i, which satisfies dfi/dt = 
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S'jHfi. Let I = (det g) 1 ^ 6 . In the isotropic case I agrees with the scale factor 
up to a multiplicative constant. Now dl/dt = HI. Hence d / '<it(/ 37 /i) = 
and the density behaves like Z -37 . From an asymptotic expression for H 
we can obtain a corresponding expression for I and hence an asymptotic 
expression for \i. For example, for the solutions whose asymptotic behaviour 
was analysed in the last section the leading order term in I is proportional 
to exp((logt)( 2 - m )/ m ). 

Consider now the case of collisionless matter with a potential satisfying 
the conditions of part a) of Theorem 4. 

Now UijU 1 ^ decays faster than any power of t and is, in particular, integrable 
in t. It is possible to argue just as in ^H] to conclude first that gij/l 2 and 
its inverse are bounded and then that gij/l 2 converges as t —>■ oo to some 
quantity gfj. It follows that 

9ij = l2 (9ij + hj) (32) 

where hij = 0(t~ n ) for any n. It is then possible to proceed as in ^H] to 
obtain asymptotics for geodesies, i.e. for the characteristics of the Vlasov 
equation. The results are strictly analogous to those in [THj except that e~ Ht 
is replaced by l~ n in the error terms in jT^j, where n is arbitrary. For a specific 
potential where the form of I can be determined these error estimates can 
be made sharper. If Vi are the spatial components of a tangent vector to a 
timelike geodesic then they converge faster than any power of t to a constant 
as t — > oo. The energy density decays like /~ 3 as in the case of dust and all 
other components of the energy-momentum tensor in an orthonormal frame 
decay faster. 



7 Conclusion 

Combining the results of this paper with previous work gives a rather detailed 
picture of the dynamics of spatially homogeneous spacetimes with accelerated 
expansion driven by a minimally coupled nonlinear scalar field with potential. 
Potentials with a positive lower bound are covered in [22], potentials which 
are positive and tend to zero at infinity slower than exponentially are treated 
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in this paper. Exponential decay with sufficiently small exponent has been 
handled in ^3], [Hj and ^H]- In the borderline case where the potential 
asymptotically approaches an exponential one but is not exactly exponential 
it would be desirable to get more information. In this paper accelerated 
expansion was demonstrated for this case for a scalar field alone in a spatially 
flat isotropic model but isotropization in more general homogeneous models, 
which is presumably true, was not. If the potential is zero somewhere then 
the expected behaviour is very different. In that case accelerated expansion 
is a transient phenomenon and the late time behaviour must be expected to 
depend in a complicated way on the Bianchi type. Cf. the discussion in the 
last section of |22j . 

The questions of global existence and the asymptotics of matter fields 
have been answered in some special cases. A treatment of more features of 
their asymptotic behaviour and of wider classes of matter fields would be 
desirable. Another direction in which the results should be generalized is to 
other mechanisms for producing accelerated expansion such as /c-essence [3]. 
See also the discussion in [23] . 

For potentials which do lead to late-time accelerated expansion an obvious 
next step is to look at inhomogeneous models. Expansions for spacetimes 
without symmetry were written down for the case of a cosmological constant 
by Starobinsky [2E| and in the case of exponential potentials in jTHJ. In [23] 
the results of j2H] were given a rigorous mathematical interpretation in terms 
of formal series and the existence of a large class of solutions (depending on 
the maximum number of free functions) was demonstrated in the vacuum 
case. Some of these results have been generalized to the case of a potential 
with a positive minimum by Bieli . In the vacuum case it was shown with 
the help of results of ^U] that all solutions evolving from small perturbations 
of de Sitter initial data have late-time asymptotic expansions of the form 
given in [2E|- Presumably a corresponding result could be obtained in any 
even spacetime dimension using the results of [3J 

For certain classes of spacetimes with high symmetry (spherical, plane 
and hyperbolic symmetry) it has been shown in and |30] that for many 
solutions of the Einstein- Vlasov system with positive cosmological constant 
there is late-time accelerated expansion and the geometry has leading-order 
asymptotics as in [2H]- It would be interesting to see to what extent these 
results can be generalized to nonlinear scalar fields with potentials which 
give accelerated expansion in the homogeneous case. The results of 31 
on the Einstein- Vlasov system coupled to a linear scalar field are a first 
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step in this direction. Note that the use of collisionless matter in these 
results is not an accident. If dust or perfect fluid with pressure were used 
then it is to be expected that there would be no global existence of classical 
solutions. Cf. the results of 21^ and [23]. The effect of nonlinear scalar 
fields on cosmological expansion in the inhomogeneous case has been studied 
heuristically and numerically in ^T]. One issue of particular interest which 
has no analogue in the homogeneous case is the formation of domain walls 
for potentials with more than one minimum. 

The way in which information about inhomogeneities in the universe is 
compared with theoretical models is via linear perturbation theory. These 
calculations are at present out of reach of rigorous mathematical theorems. 
There is not a single result proving that an inhomogeneous solution of the 
full nonlinear equations is well approximated by the expressions coming from 
linear perturbation theory. 

Cosmological solutions of the Einstein equations with accelerated expan- 
sion are a rich source of problems for mathematical relativity. It is to be 
hoped that the mathematical theory will soon have reached the point where 
it can contribute insights for the astrophysical applications of these solutions. 
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